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We consider the conditions under which solitary waves can exist in elongated clouds of Bose-Einstein 
condensed atoms. General expressions are derived for the velocity, characteristic size, and spatial 
profile of solitary waves, and the low- and high-density limits are examined. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40.Db 



Clouds of Bose-Einstein condensed atoms in elongated 
traps provide excellent conditions for investigating the 
propagation of essentially one-dimensional sound pulses 
under [jjj. In previous work, propagation of pulses in 
such traps was considered in the Thomas-Fermi approx- 
imation, and it was demonstrated that pulses propagate 
at a speed of {nUo/ra) 1 / 2 in the linear regime. Here, 
Uq = Airh 2 a sc / m is the effective two-body interaction 
matrix element, n is the particle density averaged over 
the cross section of the cloud ^,[|, m is the atomic mass, 
and a sc is the scattering length for atom-atom collisions. 
Non-linear effects were also investigated in Ref. |J and 
were found to be important for conditions of experimen- 
tal relevance. The effects of dispersion were neglected in 
this study since the length scales of interest were much 
larger than the superfluid coherence length, which sets 
the scale on which dispersive effects become important. 
It is of interest to include the effects of dispersion and 
to investigate the possible existence of solitary waves in 
these systems. 

Solitary waves have been studied in many different 
physical contexts We shall not attempt to review 

the extensive literature for the Gross-Pitaevskii equation, 
the non-linear Schrodinger equation that describes mo- 
tion of a Bose-Einstein condensate, apart from mention- 
ing the pioneering work of Zakharov and Shabat Q on 
solitary waves in bulk systems. There are a number of 
studies for clouds of Bose-Einstein condensed atoms in 
a trap. Morgan et al. considered a class of solitary 
wave solutions which are, in effect, center-of-mass oscil- 
lations of the stationary states of the time-independent 
Gross-Pitaevskii equation. 

It is possible experimentally to make very elongated 
clouds of Bose-Einstein condensed atoms in magnetic or 
optical traps, and in this paper we examine the possibil- 
ity of solitary waves propagating along the axis of such a 
cloud. The solutions we seek are the analogs of the soli- 
tary waves first observed by Scott Russell on his historic 
horse ride along the canal However, an important 
difference between canals and Bose-Einstein condensed 
clouds is that, whereas it is often a good approximation to 
assume that the depth of water in a canal is constant and 



thus the sound speed is independent of position, this is a 
poor approximation in Bose-Einstein condensed clouds, 
where the density, and hence also the sound speed, vary 
significantly in directions transverse to the axis of the 
trap. Consequently, the conditions under which one can 
find solitary wave solutions for Bose-Einstein condensed 
clouds is more restricted than for canals. We wish to con- 
sider in this paper situations where the variation of cloud 
properties along the axis of the trap may be neglected. In 
the directions perpendicular to the axis, we shall assume 
that the cloud is confined by a harmonic trapping poten- 
tial. However, to set the stage and establish notation, let 
us first consider the case of a homogeneous medium and 
then generalize our discussion to atoms in traps which 
confine particles in two directions. 

The condensate wavefunction "F satisfies the Gross- 
Pitaevskii equation, 

where V is the external potential, which may be taken to 
be zero for a homogeneous medium. We write the con- 
densate wavefunction as \& = n 1 / 2 e lc ^ > , where n is the 
particle density, which is spatially uniform in equilib- 
rium and in the absence of an external potential. To 
obtain equations resembling those of hydrodynamics, it 
is convenient to work in terms of a superfluid velocity, 
v = h'Vc^/m. We consider solutions that depend on a 
single coordinate, z. From the real and imaginary parts 
of the Gross-Pitaevskii equation one obtains two equa- 
tions, the equation of continuity 



and the equation for the phase, 

h d<t> 1 a , h 2 dW 2 

When differentiated with respect to z, the latter gives 
the generalized Euler equation, 
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where we have introduced the chemical potential, \x. For 
a dilute Bose gas, a(n) = uUq. We seek solutions to these 
equations for which the fluid velocity and particle den- 
sity propagate at a uniform velocity, u, without change of 
form, i.e., they depend on z and t only through the combi- 
nation z — ut. Thus, one may write dv/dt = —u(dv/dz) 
and dn/dt = —u{dn/dz). Far away from the solitary 
wave, the condensate is at rest and has its equilibrium 
density, uq. With these boundary conditions, the conti- 
nuity equation, Eq. (||), gives 



v = u ( 1 — 



"o 



Also, Eq. (0) can be written as 
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We integrate this equation to obtain 
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Here, we have added the integration constant —uqUq to 
impose the boundary condition n — > n at infinity. Com- 
bining Eqs. (JsJ) and (^), we obtain a differential equation 
for n. If we multiply this equation by dn 1 / 2 /dz and in- 
tegrate with respect to z, we find 
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We note that the phase has the more general form 
4>(z, t) — 4>\{z — vt) + 4>2{t), which can be seen from the 
equation dcf)/dz — mv/h. Thus, Eq. (||), can be written 
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Since the left side of this equation is a function of t and 
the right is a function of z — ut, each must be equal to 
a constant. This constant must equal —noUo in order to 
satisfy the boundary condition n — > no at infinity. Thus, 
4>2 (t) — ~ n oUot is a linear function of time. 

We see from Eq. (@) that the condition nUo — mu 2 > 
must be satisfied in order to obtain real solutions. In 
other words, the density must exceed the minimum value 



(10) 



To obtain solutions that are localized in space, n must 
lie between n m i n and no- This has a ready interpretation 
in terms of one-dimensional motion of a classical particle 



whose spatial coordinate is proportional to n 1 / 2 if z is 
regarded as the time variable. The classical potential is 
then proportional to — (n — n m i n )(n — no) 2 , and the soli- 
tary wave solutions correspond to an oscillation of the 
"particle" from n — hq to n — n m i n and back again. 
Thus, solitary waves for this problem are depressions in 
the density. In contrast, solitary waves in canals corre- 
spond to elevations of the surface of the water. We also 
note that the velocity of the wave is equal to the sound 
speed at the minimum density, n m j n . The sound speed in 
a uniform gas with density no is given by c 2 = noUo/m. 
It follows from Eq. @ that 
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Integrating Eq. (§]) we find that the profile of the wave 
is given by 



n{z) 



+ (n - n min )tanh 2 (z/C) , 



(12) 



where C = 2 1//2 £(no)[l — (n min /no)] _1 / 2 and £(n ) is the 
coherence length corresponding to the background den- 
sity no, £(no) = (87rnoa sc ) -1 / 2 . Therefore, £, which gives 
the spatial extent of the solitary wave, is on the order 
of the coherence length that corresponds to the back- 
ground density. Figure 1 shows the profile of two solitary 
waves, n(z), for Na atoms with a background density of 
n = 10 13 cm~ 3 , and for two values of n min /n , 20% 
(solid curve) and 80% (dashed curve). 
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FIG. 1. Density profiles of solitary waves as function of 
space, as given by Eq. ([[2]). Here the background density 
of atoms is no = 10 13 cm -3 and therefore for Na atoms 
£(no) ~ 1.2 ;Um. For the solid curve, the ratio n m i n /no = 
20% and £ w 1.9 /im; for the dashed curve, n m i n /no = 80% 
and £ ?s 3.8 /im. 
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For typical densities in the MIT experiment (n « 10 14 
cm -3 at the center of the cloud), £ is ~ 0.2 — 0.4 /jm 
with a sc = 27.5 A for Na atoms |^|. However, in or- 
der to be able to observe solitary waves it is desirable 
to look at lower densities than these, so that the coher- 
ence length becomes sufficiently large that structures can 
be resolved by optical means. For non-zero n m i n the so- 
lutions ([l2]) correspond to what are termed "gray soli- 
tons" in Ref. IK}. Finally, we point out that, if u = 

(^min 

= 0), EqT^12j) gives the well-known kink solution 
= n, 1 / 2 tanh(z/e) JlT|, which is sometimes referred to 



as a "dark soliton" (lO[0. 

We turn now to the more realistic problem of atoms in 
a trapping potential which is harmonic in the transverse 
direction and for which there is no restoring force along 
the axis of the trap, which we take to be the z axis. We 
assume that the transverse dimension of the cloud is so 
small that the time scale for adjustment of the transverse 
profile of the particle density to the equilibrium form ap- 
propriate for the instantaneous number of particles per 
unit length is small compared with the time for the pulse 
to pass a given point. Later, we shall investigate what 
this condition means quantitatively. The problem be- 
comes one-dimensional, and the solitary pulse may be 
specified in terms of a local velocity, v(z), and a local 
density of particles per unit length, a(z), 



a(z) = J dxdy\ i $(x,y,z)\ i 



(13) 



Here, x and y are coordinates perpendicular to the axis 
of the trap. With this assumption, the wavefunction may 
be written in the form 



f(r,i) = f(z,t)g(x,y,a) 



(14) 



where g is the equilibrium wavefunction for the transverse 
motion; g(x, y, a) depends on time implicitly through the 
time dependence of a. We choose g to be normalized so 
that J \g\ 2 dxdy = 1, and therefore from Eqs. (O) and 
@,J/| 2 = a. 

It is convenient to derive the equations for / and g 
from a variational principle. The Ginzburg-Pitaevskii 
equation may be derived by requiring stationarity of the 
action 
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Using Eq. (|l4| ) for we can write Eq. (|T^ ) as 
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where V 2 ^ = d 2 /dx 2 + d 2 /dy 2 . Terms containing 
(df /dz)(dg/dz) vanish because of the normalization con- 
dition on g. Minimizing Eq. ( [l6| ) with respect to g* , we 
find that g obeys the equation 
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where [i is the chemical potential. We neglect the last 
term on the left side of Eq. (|l7|) since, as we show below, 
the characteristic length of pulses is sufficiently long that 
it is negligible in all cases of interest. We now minimize 
Eq. (|l6l ) with respect to /* and find 
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+U a (J \g\ 4 dxdy) |/| 2 / + (J \g\ 2 Vdxdy) f . (18) 

Following the same procedure as for the homogeneous 
medium, we assume that / = <r 1 / 2 e 1 ^, with the velocity 
field associated with / being v = (h/m)dcj)/dz. Again, 
we obtain hydrodynamic equations for a and v, which 
are the same as Eqs. (||) and (||) but with n replaced by 
a and /x given in the present case by 
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The first term on the right side of Eq. (|19|) is the kinetic 
energy in the transverse direction, the second is the po- 
tential energy due to the confining potential, and the 
third is the energy due to interactions between atoms. 
Assuming that a and v are functions of z — ut only with 
u constant, we find that 
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where e (a) — (i(a')da' is the energy per unit length. 
In Eq. (|20|), we have imposed the boundary conditions 
a = (Jq and v = far away from the disturbance. 

There is no simple expression for the energy per unit 
length for general values of a. To obtain analytical re- 
sults, we explore some limiting cases for the experimen- 
tally relevant situation where the confining potential in 
the transverse directions is harmonic and rotationally in- 
variant, and given by V = mLo 2 ± {x 2 + y 2 )/2. In the 
low-density regime, the interaction energy can be treated 
perturbatively and the problem reduces essentially to the 
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one-dimensional case treated above. As we shall show be- 
low, in this limit e(cr) oc er(l + a sc a). In the high-density 
limit, the Thomas-Fermi approximation applies jl3), and 
e(cr) varies as cr 3 / 2 . 

Before examining the two cases separately, we estimate 
the density, cr c , at which the cross-over between the two 
limits occurs. This corresponds to the condition that 
the interaction energy per particle is on the order of the 
oscillator energy in the transverse direction. The inter- 
action energy per particle is nUo. If A is the cross sec- 
tion of the cloud and R± is the corresponding radius, 
n = a I A. Therefore, the condition determining a c is, 
Uq(t c /A ~ huj±. Denoting the characteristic length scale 
for the ground state of a particle in the transverse con- 
fining potential by a± = (h/mui±) 1 ^ 2 and assuming that 



R± = a± , the cross section of the cloud is A 



thus 



one finds a c ~ a" 1 , which gives a c 4 x 10 6 cm -1 for Na 
atoms. Alternatively, we can determine a c by equating 
the interaction energy per particle and the kinetic energy 
of the atoms due to their confinement in the transverse 
direction, a c U /A ~ fi 2 /(mA), which gives the same re- 
sult for a c . This expression for a c is equivalent to the 
condition that the coherence length, £, be comparable to 
the transverse dimension of the cloud, R± . 

We now examine the low-density limit, a <C a c in 
which the interaction energy can be treated perturba- 
tively. To find the differential equation that g satisfies 
in this limit, we neglect the interaction-energy term (the 
third term on the left side) in Eq. (JTt|) . The last term 
on the left side is ~ h 2 / (2m( 2 )(5f /fy , where Sf is the 
deviation of /, which is much less than h 2 /(2ma 2 L ) be- 
cause of the condition £ 3> a±. Therefore, this term is 
negligible, and g thus satisfies the Schrodinger equation 
in a harmonic potential, 
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g + Vg = fig 



(21) 



which means that \g\ oc exp [— (x 2 + y 2 )/2a 2 L \. From 
Eq. jl9|), we find that fx(a) = hu>±(l + 2a sc a). The 
first term in /i is the ground-state energy of the har- 
monic oscillator and the second is the interaction energy. 
Thus, e(cr) = %uj±_u{\ + o sc a) in the low-density limit, 
and Eq. (EG) takes the form 
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by analogy with Eq. (g). The results for the profile 
and the width of the cloud are analogous to the one- 
dimensional case, u 2 = a mul Uo/(2mA) 1 where a m - m is 
the minimum of a associated with the solitary wave and 



OT C = 2 1 /2£( no /2)[l - (<r miD /a )}-^ 2 with n = (To/A, 
From Eq. ( |22| ) , we see that the density which determines 
u 2 is er m i n /2A and the one that determines £ is ao/2A. 
The factor of 1/2 in these results compared with the anal- 
ogous results for the homogeneous case is due to the av- 
erage of the equilibrium density of atoms across the trap, 
n(x,y), given by J n 2 (x , y) dxdy / J n(x , y) dxdy . Its ori- 
gin is thus completely different from that of the factor of 
1/2 that occurs in the Thomas-Fermi limit Q. The latter 
is J n(x,y) dxdy /[n(0,0) J dxdy}. It is amusing that the 
sound speed expressed in terms of the density of parti- 
cles on the axis of the trap is given by precisely the same 
expression, c = [n(0, 0)Uo/(2m)] 1 ^ 2 , in both the high- 
and low-density limits for harmonic transverse trapping 
potentials. 

We argued above that the motion would be quasi-one- 
dimensional if the time scale for adjustments in the trans- 
verse structure of the cloud is short compared with the 
time for passage of the pulse, and we now check the con- 
sistency of this assumption. In the low-density regime, 
a <C <7 C , the characteristic time for adjustment of the 
profile is on the order of wj , while the time scale for 
passage of the pulse is at least of order £/c. The ratio of 
these scales is nUo/hLO±, which is much less than unity in 
this regime, and therefore our assumption is consistent. 

Another case that can be solved analytically is that 
of small-amplitude solitary waves for arbitrary densities. 
We expand e(<r) to order (a — ctq) 3 and write Eq. (|20| ) as 
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where the prime denotes differentiation with respect to a. 
Clearly, two roots of the right side of the above equation 
are equal to ctq. The remaining root gives an expression 
for the velocity u as a function of Ctq an( i Cmim 



Oo) + 



(ffmin — CT ) 
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Because e (cr) depends on a, the velocity of propagation 
of the wave is not generally equal to the bulk velocity 
at the minimum density, except in the low-density limit, 
when e (a) is independent of a. Integrating Eq. (|25|), we 
can determine the profile of the solitary wave, a{z). For 
small values of a — ctq, Eq. (p5|) takes the simple form 



da 1 ' 2 (a - cr min )V2 



dz 
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(27) 



cr(z) 



+ (ctq - (Tmin) tanh 2 (z/C) . 



For this problem, 

C = 2e(n )[l-(cr min /cr )]- 1 /2 



(23) 



(24) 



where the length I is given by I 2 — h 2 /(mao[e (ctq) + 
"of (cr )/3]). Again assuming that a is close to ctq, we 
find that 
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!)tanh 2 (z/C) 
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with C = l[l - OwAro)]- 172 - 

We now consider the calculation of e(a) in the high- 
density limit where the kinetic- energy term is negli gibl e 
fL3"f . We also neglect the last term on the left of Eq. (17). 
The consistency of this assumption will be checked below. 
With these approximations, we write the Thomas-Fermi 
equation for g, 



Vg + U \fg\ 2 g = M 



(29) 



Thus, \g\ 2 is a parabola, \g\ 2 oc [1 - {x 2 + y 2 )/R 2 ± }- This 
form for \g\ implies that the potential energy is equal 
to muj 2 L R 2 _/Q and that the interaction energy is equal 



to 4U a/(3irR 2 ± ). From Eq. ( |19| ) we see that /i(a) is the 
sum of these two terms in this limit. To find the explicit 
dependence of [i on c, we calculate R± |j| now. The 
density of atoms, n(x,y,z), has the functional form of 

M 2 , 



n{x, y, z) = n(0, 0, z) ( 1 



(30) 



where n(0, 0, z) is the density on the axis of the trap. 
Thus, the number of particles per unit length is given by 

<j{z) = J n(x, y, z) dxdy = — n(0, 0, z)~kR\ . (31) 

The density on the axis of the trap is given by 



u(0,0, z)U = -rmj\R\ 



(32) 



Thus, we find that R 2 ± = Aa\ (aa^) 1 / 2 . Equation (|l|) 
then implies that yu.(tr) = 2hLUj_(aa ec ) 1 ^ 2 and therefore 
e(a) = {4/3)hu±(a sc a 3 ) 1/2 . Going back to the width, f, 
of the solitary waves, we find that C = (12/5) 1 / 2 ^(n )[l- 
(W^o)]" 1 ' 2 where n = cr / {tt R 2 ^) . 

In considering the limits of validity of our calculation 
in the high-density regime, we observe that the time scale 
for adjustment of the profile of the pulse is ~ R±/c and 
that the time for passage of the pulse is ~ £/c. If the 
motion is to be essentially one dimensional, R± must be 
much smaller than £. In the Thomas-Fermi approxima- 
tion, R± is larger than a± due to the repulsive interac- 
tions between the atoms. This can be seen from the for- 
mula for the radius R± — 2a± (erase) 1 / 4 derived above. To 
satisfy the condition £ 3> R± , the quantity aa sc must be 
large compared with (3/10) 2 (croa sc ) _1 [1 — (oWn/co)] -2 - 
Thus, if da = oq — cr m in is the amplitude of the distur- 
bance, we see that Sa/aa <C 3/(10<7oa sc )- This indicates 
that the amplitude of the solitary wave, Sa, must be ex- 
tremely small if the one-dimensional approximation is to 
be valid when a 3> a c , and therefore our small- amplitude 
treatment given above is applicable. For the experimen- 
tal conditions of Ref . [ lj , the number of particles in the 
trap is N ~ 5 x 10 6 , and the length of the trap in the axial 
direction L is ~ 450 /im, implying that <7o ~ 10 8 cm" . 
Validity of the one-dimensional approximation therefore 



requires that 8a/ ao <C 0.01. Finally, we check the consis- 
tency of the assumption made in deriving the Thomas- 
Fermi equation for g. The last term on the left of Eq. ( |l7| ) 
has an upper bound of Ti 2 / (2m£ 2 ). This term is approxi- 
mately equal to h 2 /(2m£ 2 )(5a/a ) < h 2 /(2m£ 2 ) ~ nU . 
Therefore, this term is indeed negligible. 

Our calculations indicate that the most favorable con- 
ditions for observing solitary waves in trapped Bose- 
condensed gases occur when the density is sufficiently 
low to permit a perturbative treatment of particle inter- 
actions. In this case, one-dimensional behavior persists 
even for large- amplitude solitary waves. Low-density 
systems have the further advantage that the coherence 
length is correspondingly large. Since the coherence 
length determines the size of solitary waves, this sim- 
plifies resolution of these structures. However, a lower 
density of particles makes the detection of these effects 
more difficult because of the lower signal. 

We now wish to estimate the experimental conditions 
required for observation of solitary waves in the low- 
density regime. The spatial resolution of current experi- 
ments using direct imaging methods is ~ 4 /im ]l4} |. As 
a theoretical estimate for the width of the solitary wave, 
we take the full width of the dip in the density profile 
at half the maximum depression, that is the distance be- 
tween points where a = (ao + er m i n )/2. This is given by 
2Ctanir 1 (l/V2) » 3.5£(n )[l - (a mi „/a )]- 1/2 , which 
must exceed the experimental resolution if solitary waves 
are to be observable. This leads to the condition 



£(n )>[l-(<W(7o)] 1/2 /im 



(33) 



or £(no)x^l /jm for solitary waves with a mul /ao *C 1. 
For the coherence length to exceed 1 /urn, the density per 
unit volume must be less than n b s ~ 10 13 cm~ 3 . The 
corresponding number of particles per unit length a is 
n Q bs^4 w riohsTra 2 ^. For a transverse trapping frequency 
of 240 Hz, which is a typical value of the MIT trap |],[01 , 
the oscillator length, a±, in this direction is ~ 1.4 //m, 
implying that the number of particles per unit length 
must be less than 8 x 10 5 cm -1 for structures to be ob- 
servable. This value corresponds to a c /5, and therefore 
the coupling would be weak. Let us now calculate the 
number of particles which should be used in the MIT trap 
for a to have the value 8 x 10 5 cm -1 . The kinetic energy 
along the axial direction is typically of order h 2 /2mL 2 , 
and the interaction energy is ~ aUo/(icaj_). The ra- 
tio of these quantities is ~ (a^/L) 2 (8aa sc )~ 1 . Thus for 
L ^> ax/(8cra sc ) 1 / 2 , the kinetic energy along the z direc- 
tion is negligible, and the Thomas-Fermi approximation 
]l3t may be used to determine the structure along the 
axis of the trap. We shall assume this to be the case and 
subsequently will check the consistency of this assump- 
tion. Under these conditions, the Thomas-Fermi approx- 
imation may be used to calculate the z-dependence of a, 
even though it may not be used to calculate the struc- 
ture of the cloud in the transverse directions. In the 
presence of a confining potential in the z direction, the 
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Thomas-Fermi condition is that the sum of the chemi- 
cal potential and the z-dependent part of the trapping 
potential should be constant. Since the chemical poten- 
tial is fj, = hu)±(l + 2aa sc ) in the low-density limit, this 
condition is 



Tiu!±[l + 2a(z)a sc ] + ^mu> 2 z 2 = constant, 



(34) 



where lu z is the frequency of the trapping potential al ong 
the z axis. We equate the value of the left side of Eq. ( |34| ) 
at z = to the value at z — Z, where Z is the distance 
from the center of the cloud to its edge, Z = L/2. Since 
a vanishes at the edges of the cloud, we get 



2htu±a(z = 0)a s . 



1 



^ 2 ryl 



(35) 



where a(z = 0) is the number of particles per unit length 
at the center of the cloud. Solving the above equation 
for Z, we find 



= 2^Hz = Kc ] 1 / 2 
a± 

= 2^{n ohs a sc ) 1/2 a 2 z , 



(36) 



where a z — (Ujmuj z ) 1 l' 1 is the oscillator length in the 
axial direction. From Eqs. (34) and (|36| ) we see that the 
number of particles per unit length can be written as 



a(z) =a(z = 0)(l 



z 

'z 2 



(37) 



The total number of particles is N = f_ z a(z)dz; thus 
from Eqs. @ and (§7|) we find 



N = -a(z = 0)Z 



8tt 3 / 2 



n% 2 s (a z a ± ) 2 aU 2 ■ (38) 



For the MIT trap, the frequency in the axial direction 
is ~ 7.9 Hz |Tj| and therefore a z ~ 7.5 fim; the total 
number of particles, N, should be equal to ~ 2 x 10 3 
if a(z = 0) is to be 8 x 10 5 cm" 1 . The total width of 
the trap in the z direction would then be L as 38 /im. 
Returning to our assumption, we see that for a = er c /5 
(the highest value of a for which these structures are 
observable), the condition L ^> a±/(8aa sc ) 1 ^ 2 implies 
that L must be much larger than 0.8a±. This is indeed 
true and the Thomas-Fermi approximation is valid along 
the z axis of the cloud. 

The development of traps with stronger transverse con- 
finement, such as the optical dipole trap of the MIT 
group JlTjf , promises to facilitate experiments on solitary 
waves since, for a given value of the dimensionless cou- 
pling <7<Zsc, the corresponding densities will be higher. 
However, the disadvantage is that the higher density will 
result in shorter length scales for structures. 

In this paper we have discussed the propagation of soli- 
tary waves in elongated traps and have estimated charac- 
teristic sizes of these structures. Clearly, many questions 



remain. These include the stability of these structures, 
the role they play in other phenomena (such as dissipa- 
tion) , and the question of how pulses propagate when the 
motion is not quasi-one-dimensional. 

Helpful discussions with H. Fogedby, W. Ketterle, H.- 
J. Miesner, M. Saffman and H. Smith are gratefully ac- 
knowledged. G.M.K. would like to thank the Foundation 
of Research and Technology, Hellas (FORTH) for its hos- 
pitality. 
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